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Summary. — We present a summary of the measurements of the top quark 
production and decay properties performed by the DO experiment at the Fermi- 
lab Tevatron proton-antiproton collider at ^/s = 1.96 TeV using 1 fb~^ of data. 
We discuss the first simultaneous measurement of the ratio of branching fractions, 
R = B{t -^ Wb)/B{t — > Wq), with q being a d, s, or b quark, and the top quark 
pair production cross section au, the first model-independent measurement of the 
helicity of W bosons produced in the top quark decays, and the first measurement 
of the integrated forward-backward charge asymmetry in the production of ii pairs. 

PACS 13 . 85 . Rm - Limits on production of particles. 

PACS 13.85.Qk - Inclusive production with identified leptons, photons, or other 

nonhadronic particles. 

PACS 14. 65. Ha - Top quarks. 

PACS 14. 70. Ha - W bosons. 
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1. — Introduction 

The large sample of top-antitop quark pairs produced by the Fermilab Tevatron col- 
lider allows not only to perform precision measurements of such fundamental top quark 
characteristics as its production cross section and mass, but also to study a broad variety 
of top quark production and decay properties to address the question whether the top 
quark is indeed the particle predicted by the standard model (SM). Any deviation from 
the SM prediction would be an indication of the new physics. In this paper, we study the 
properties of top quark decay by measuring the ratio of top quark branching fractions, 
R = B{t ^ Wb)/B(t -^ Wq) [1], and the W boson hehcity [2]. We also report on a 
measurement of forward-backward (/6) charge asymmetry in tt production [3]. 

2. Selection of ti candidates 

To study top quark properties we mainly use the top quark pair decay channel where 
one W boson from a top quark decays into two quarks, and the other one into an electron 
or muon and a neutrino, referred to as the lepton plus jets (£-|-jets) channel. We select 
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a data sample enriched in tt events by requiring a lepton and a jet at the trigger level. 
We further require at least four jets with transverse momentum px > 20 GeV and 
pseudorapidity \ri\ < 2.5, one isolated electron (muon) with px > 20 GeV and \ri\ < 1.1 
(I77I < 2.0), no other isolated lepton with pT > 15 GeV, and missing transverse energy 
]^T > 20 GcV (e+jets) or ]^t > 25 GeV (yu+jets). The leading jet pr is required to 
exceed 40 GcV. 

The dominant background in the selected sample is the production of W bosons in 
association with heavy and light flavor jets (W^+jets). Smaller contributions arise from 
Z+jets, diboson and single top quark production. We model W/Z+jets production 
with the ALPGEN [4] Monte Carlo (MC) generator for the matrix element calculation 
and PYTHIA [5] for parton showering and hadronization. We use pythia and single- 
top [6] event generators to model diboson processes and single top quark production, 
respectively. The sample also includes contribution from multijct events in which a jet 
is misidentified as an electron (e+jets channel) or in which a muon originating from ei- 
ther a semileptonic decay of a heavy quark or an in-flight pion or kaon decay in a light 
flavor jet appears isolated (/x-|-jets channel). We determine multijet contribution from 
data. Details of the object identification, selection and background evaluation in £-|-jets 
channel can be found elsewhere [7] . 

To improve the purity of the samples used for top quark property measurements 
and/or discriminate ti signal from background we use information about the presence 
of b jets in the event. We identify b jets using the algorithm which combines variables 
that characterize the presence and properties of secondary vertices and tracks with high 
impact parameter inside the jet into a neural network [8]. 

We also utilize kinematic information to separate ti events in the selected samples 
from the background. We consider kinematic variables which discriminate between the 
ti signal and background and combine them into a discriminant function. The selected 
variables are required to be well described by the background model. Neglecting the 
correlations between the input variables, the discriminant function can be approximated 
by the expression [7]: 

p^ Il^S^ix^)/b,{xi) 



1 + 'n.iSi{Xi)/b^{Xi) 

where Si(xi) and bi{xi) are the normalized distributions of variable i for signal and 
background, respectively. As constructed, the discriminant function peaks near zero for 
the background, and near unity for the signal. We use the discriminant distribution cither 
to determine the fraction of ti events in a sample via a Poisson maximum-likelihood fit 
to data of a sum of signal and background discriminant distributions [1, 3], or to improve 
the purity of the samples by selecting events with high values of V [2] . 

3. Simultaneous measurement of the ratio of branching fractions R — 3(t — > 

Wb)/B(t -^ Wq) and (7^^. 

Within the SM top quarks decay almost exclusively to a T4^ boson and a b quark. 
This prediction is based on the requirements that there are three fermion families and 
the CKM matrix [9] is unitary. If these assumptions are relaxed, CKM matrix clement 
\Vtb\ is essentially unconstrained, \Vtq\ elements can significantly deviate from their SM 
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values, and the ratio R of the top quark branching fractions, which can be expressed as 



R 



B{t -^ Wb) 
B{t -^ Wq) 



Vtb 



Vtb 



Vu 



Vtd 



can be different from unity. A precise measurement of R aUows to perform direct mea- 
surement, free of assumptions about the number of quark famihes or the unitarity of the 
CKM matrix, of the | Vtq \ elements via the combination with measurements of the single 
top quark production rates in s and t channels [10]. 

The probability for a tt event to have certain number of identified h jets depends on 
the jet flavor, and therefore it depends on R. Fig 1(a) shows the fraction of events with 
0, 1 and > 2 tagged jets as a function of R for tt events with > 4 jets. If R is close to 
zero, i.e., top quarks decay predominantly to W bosons and light quarks, the probability 
to have a tt event with two b tags is negligible while the probability to have no tags is 
close to 90%. Fig 1(b) shows the comparison of the number of data events with 0, 1 or 
>2b tags to the prediction for the sum of background and tt signal with R = 0, R = 0.5 
and i? = 1 . 
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Fig. 1. - (a) Probability of ti events to have 0, 1 and > 2 6 tags as a function of R for events 
with > 4 jets; (b) predicted sum of tt signal and background for R — 0, 0.5 and 1 compared to 
data in the 0, 1 and > 2 6 tag samples for events with > 4 jets. 



We measure the ratio R using f+jets events with three or more jets. We split this 
sample into subsamples according to lepton flavor (e or /i), jet multiplicity (3 or > 4 
jets) and number of identified 6-jets (0, 1 or > 2), and extract R and the ti production 
cross section att simultaneously from the fit of the predicted numbers of ti signal and 
background events to the observed ones with 0, 1 and > 2 6 tags. For the events 
with > 4 jets and b tags, we include the shape of a topological discriminant in the 
fit. We build the discriminant function using simulated PF+jets and ti events. For the 
latter we obtain a distribution for each of the three decay modes of ti: ti -^ W~^bW~b, 
ti -^ W^bW'qi (or ti -^ W~^qiW~b) and ti -^ W~^qiW~qi, where qi denotes a hght 
down-type (d or s) quark. The systematic uncertainties are incorporated in the fit using 
nuisance parameters [11], each represented by a Gaussian term in the likelihood fit. The 
result of the fit is: 



i? - 97+009 
n — 'J-Ji -0.08 



(stat + syst) and 
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att = 8.18Ij;:^^ (stat + syst) ± 0.50 (lumi) pb 



for a top quark mass of 175 GeV. Fig 2 compares the distribution of the data to the sum 
of predicted background and measured signal for events with exactly three and > 4 jets. 
We extract a limit on R and |Vt6| following the Feldman-Cousins procedure [12] which 
yields R > 0.79 at 95% C.L. and \Vtb\ > 0.89 at 95% C.L. The latter hmit is derived 
assuming a unitary CKM matrix with three fermion generations. We also determine a 
model-independent limit of the ratio of \Vtb\^ to the off-diagonal matrix elements to be 
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Fig. 2. - Predicted and observed number of events for the measured R and au in the 0, 1 and 
> 2 6 tag samples for events with (a) three jets, (b) > 4 jets, and (c) predicted and observed 
discriminant distribution in the b tag sample with > 4 jets 



4. Measurement of the W helicity in top quark decays 

In the SM the Wtb top quark decay vertex has a V — A structure which defines the 
fractions of W bosons produced in each polarization state to be 0.697 ± 0.012 [13] and 
1.6 X 10""' [14] for the longitudinal (/g) and right-handed (/+) fractions, respectively, for 
a top mass of 172.5 GeV. Any significant deviation from these values will be a signal of 
new physics. 

We measure the W helicity by studying the distribution of the angle 9*, defined in 
the W boson rest frame, between the momentum of down-type fermion (charged lepton 
or d or s quark) from the W boson decay and the top quark. The distributions of 
cos 6** expected for the pure right-handed, longitudinal and left-handed VF's are shown 
in fig 3(a). 

The data sample used in this analysis includes -^-fjets events and dilepton events 
where both W bosons from top quarks decay into leptons (e or /i). We select such 
events by requiring two charged leptons with opposite charge and pt > 15 GeV, and 
at least two jets with pT > 20 GeV. To increase signal purity we built a topological 
discriminant in each of the five channels considered (ee, e/i, pp, //-f jets. e-|-jets) using 
the kinematic variables listed in tab 1. Variables NNi,-^ and (NNf,), the output value of 
the neural network of the leading jet and the mean of the neural network outputs of the 
two leading jets, respectively, include information about jet flavors into the discriminant. 
Thus, instead of requiring a b tag in the event corresponding to a cut on NNb, we use 
the full NNb distribution and avoid a loss of efficiency due to b tag requirement. An 
example of the discriminant distribution in the e-t-jets channel is shown in fig 3(b). In 



MEASUREMENT OF TOP QUARK PROPERTIES AT DO 5 

Table 1. - Summary of the m,ultivariate selection and number of selected events for each of the 
ti final states in W helictty measurement. The uncertainties are statistical only, except for the 
background estimates in the ee and fifi channels, m which systematic uncertainties arising from 
imperfections in the MC model of the data are included. 
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each channel we select a cut on V to achieve the best expected precision for W helicity. 
Tab 1 lists the cut values and the composition of each sample after the cut. 
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Fig. 3. - (a) Distribution of the cos^* for pure left-handed (green curve), right-handed (red) 
and longitudinal (black) W boson helicity; (b) distribution of discriminant for signal (open 
histogram) and background (green) in the e±jets channel. 



To calculate cos 0* we reconstruct events in £+jets channel using a kinematic fitter. 
The fitter varies the four-momenta of the objects in the event within their resolutions 
and minimizes a x^ under the following constraints: i) both W boson masses are required 
to be exactly 80.4 GeV; ii) the masses of two reconstructed top quarks are required to be 
exactly 172.5 GeV. Among possible 12 solutions we choose the one with the best x^ from 
the kinematic fit and with the maximim probability based on NNi, values of all four jets. 
In the dilepton channel, the kinematics are underconstrained. We reconstruct dilepton 
events assuming a top quark mass of 172.5 GeV. Each event provides two entries to the 
cos 9* distribution. Fig 4 shows the reconstructed cos 0* distributions for (a) leptonic and 
(b) hadronic W boson decays in £±jets channel, and (c) for dilepton channel. Although 
hadronic W decays do not allow to discriminate between left-handed and right-handed 
W bosons, they do add information on the longitudinal W boson fraction. Including the 
hadronic decays in the measurement improved the sensitivity by 20%. 

We extract /o and /+ from the binned Poisson likelihood fit of the sum of background 
and right-handed, left-handed, and longitudinal helicity templates for signal to data. The 
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measured values of /o and /+ are: 

/o = 0.425 ± 0.166(stat) ± 0.102(syst) 
/+ = 0.119 ± 0.090(stat) ± 0.053(syst) 

with a correlation factor of -0.83. Fig 5 shows the result of the fit along with the 68% 
and 95% C.L. contours from the fit, including all uncertainties. The measured values are 
consistent with the SM, as there is a 30% chance to observe a larger discrepancy given 
the statistical and systematic uncertainties of the measurement. If /+ is fixed to the SM 
value, we find /o = 0.619 ± 0.090(stat) ± 0.052(syst) and if /o is fixed to the SM value, 
/+ yields /+ = -0.002 ± 0.047(stat) ± 0.047(syst). 
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Fig. 4. - Comparison of the cos 6* distributions in data (points with error bars) with the result 
of the fit (solid open histograms) and the standard model expectation (dashed open histograms) 
for (a) leptonic W boson decays and (b) hadronic W boson decays in ^-fjets evens, and (c) in 
dilepton events. 



5. Measurement of the forward-backward charge asymmetry in tt production 

Leading order calculations of ti production in pp collisions predict that the kinematic 
distributions of a top quark pair are charge symmetric. But this symmetry is acciden- 
tal. At the higher orders top pair production is expected to be charge asymmetric. In 
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Fig. 5. - Result of the model-independent W boson helicity fit. The ellipses show the 68% and 
95% C.L. contours, the triangle corresponds to the physically allowed region where /o + /+ < L 
The star denotes the SM prediction. 
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particular, at NLO the predictions range from 5% to 10% [15]. NNLO calculations [16] 
are not complete but deviate significantly from the NLO predictions. But in both cases 
expected asymmetry is small. If new physics is present in tt production, the observed 
asymmetry can be much larger than that predicted by the SM, thus making the charge 
asymmetry measurement a probe for new physics. 

A charge asymmetry in ti production can be observed as a forward-backward inte- 
grated asymmetry defined as Aft = {Nf — Ni,)/{Nf + Nt), where Nf{Ni,) is the number 
of events with a positive (negative) Ay and Ay = yt — yt is the signed difference between 
top and anti-top quark rapidities. The studies using events simulated with MC@NLO 
event generator [17] showed that integrated asymmetry strongly depends on the region 
of phase space under study. In particular, we observe large variations of asymmetry from 
-1-8% to -3% as a function of the fourth highest particle jet px- It also shows a strong 
dependence on the number of jets in the event. 

For this measurement we select £-|-jets events with at least four jets, one of which is 
required to be identified as a 6 jet by the neural network algorithm. We reconstruct the 
kinematics of the top quark pairs using the kinematic fitter described above with the 
difference that top quark mass is assumed to be 170 GeV. We use the 6-tagged and the 
three remaining highest pt jets in the fit, thus reducing the number of possible parton-jet 
assignments considered in the fit. We select the solution with the lowest x^- 

We build a discriminant to separate ti signal from the background using the following 
variables: i) the px of the leading 6-tagged jet, ii) the x^ of the kinematic fit. Hi) the 
invariant mass of two jets assigned to W boson by the fit, and iv) kJp"^ — p™"'-R™™, 
where i?™*" is the smallest angular distance between any two jets used in kinematic fit, 
and p™'" is the smaller of these two jets' transverse momenta. 

We extract the ii content and Afh simultaneously by performing a maximum-likelihood 
fit of the sum of backgrounds and backward and forward signal to the discriminant dis- 
tribution in data and to the sign of the reconstucted Ay. Forward and backward signal 
have the same discriminant distribution but different sign of the reconstructed Ay. The 
fitted sample composition compared to the discriminant distribution in data for events 
with > 4 jets is shown in fig 6. We measure asymmetry for events with > 4, = 4 and > 5 
jets. We find Afi, = (12 ± 8)%, Afb = (19 ± 9)%, and Afb = (-16li^)%, respectively, 
consistent with the NLO predictions from MC@NLO. 
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Fig. 6. - Distribution of the discriminant for data and a sum of tt signal and background for 
events reconstructed as forward (Aj/reco > 0) (a) and backward (Ay,.eco < 0) (b). 



To compare the measured integrated asymmetry with theory we either have to correct 
it for the reconstruction and acceptance effects or provide a prescription that would allow 
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such comparison. We have chosen the latter option since correction can reduce sensitivity 
to the effects beyond the SM. We define the asymmetry at the particle level as 

g{\Ay\) - g{-\Ay\) 
^^^^'^^'^ ^ 9i\Ay\) + gHAy\) ' 

where g is the probability density for generated |Aj/| within the acceptance. Since the 
measured asymmetry can be diluted due to misreconstruction of the event kinematics or 
misidentification of the Icpton charge the predicted reconstructed asymmetry is given by 

/•OO 

A^'" - J^ Af,{Ay)D{AyMAy) + g{-Ay)]dAy 

where D is the dilution function. It is defined as I? = 2P — 1, where P is the probability 
to reconstruct the correct sign of Ay. This function is determined using simulated tt 
events passed through the full reconstruction chain. Dilution represents the fraction of 
asymmetry visible in the detector and is parameterized as 

Di\Ay\) = colnil + ci\Ay\ + C2\Ayf) , 

where the parameters depend on the number of jets in the event. The dependence of the 
dilution function on \Ay\ introduces model dependence into any correction of Afi, from 
the observed to the particle- level asymmetry as they depend on the model's Ajb(|Aj/|) 
shape. An example of the dilution function is shown in fig 7(a) and the parameters can 
be found in [3]. 

A contribution from the new physics can affect the ti forward-backward charge asym- 
metry. Here we consider the sensitivity of the measurement to tt production via heavy 
neutral gauge bosons {Z') decaying predominantly to quarks. We study the scenario 
where the coupling between the Z' boson and the quarks is proportional to that of Z 
boson and the quarks, and the interference effects with the SM tt production are negli- 
gible. Unlike the direct searches of Z' [18] the asymmetry measurement is sensitive to 
the production of both narrow and wide resonances. We simulated ti events produced 
via narrow Z' resonance as in [18] using pythia and found large positive asymmetries 
ranging from 13% to 35% depending on the resonance mass. We predict the distribution 
of Afb as a function of the fraction of ti events produced via Z' for each Z' mass and set 
limits following the Feldman-Cousins procedure [12] (see fig 7(b)). 

6. Summary 

In summary, we have presented the measurements of the production and decay prop- 
erties of the top quark. All measurements are consistent with the SM. The large data 
set collected at the Tevatron will allow more precise measurements of the top quark 
properties in the near future and facilitate searches for new physics in the top quark 
sector. 
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Fig. 7. - (a) The geometric dilution and its uncertainty band as a function of generated \Ay\ 
for the standard model tt production and > 4 jets; (b) 95% C.L. limits on the fraction of ti 
produced via a Z' resonance as a function of the Z' mass. Expected (observed) limits are shown 
by the dashed (solid) curve. Shaded bands correspond to one and two standard deviations from 
the expected limit. The excluded region is hatched. 



